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Abstract 

In the present paper, Unruh-DeWitt detectors are used in order to investigate the 
issue of temperature associated with a spherically symmetric dynamical space-times. 
Firstly, we review the semi-classical tunneling method, then we introduce the Unruh- 
DeWitt detector approach. We show that for the generic static black hole case and 
the FRW de Sitter case, making use of peculiar Kodama trajectories, semiclassical and 
quantum field theoretic techniques give the same standard and well known thermal 
interpretation, with an associated temperature, corrected by appropriate Tolman fac- 
tors. For a FRW space-time interpolating de Sitter space with the Einstein-de Sitter 
universe (that is a more realistic situation in the frame of ACDM cosmologies), we 
show that the detector response splits into a de Sitter contribution plus a fluctuat- 
ing term containing no trace of Boltzmann-like factors, but rather describing the way 
thermal equilibrium is reached in the late time limit. As a consequence, and unlike the 
case of black holes, the identification of the dynamical surface gravity of a cosmologi- 
cal trapping horizon as an effective temperature parameter seems lost, at least for our 
co-moving simplified detectors. The possibility remains that a detector performing a 
proper motion along a Kodama trajectory may register something more, in which case 
the horizon surface gravity would be associated more likely to vacuum correlations 
than to particle creation. 

PACS: 04.70.-S, 04.70.Dy 



1 Introduction 

It is well known that Hawking radiation [T] is considered one of the most impor- 
tant predictions of quantum field theory in curved space-time. Several deriva- 
tions of this effect have been proposed and recently the search for "exper- 
imental" verification making use of analogue models has been pursued by many 
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investigators (see for example [ZIIH])- One of the most beautiful achievements is 
represented by the interpretation of the surface gravity k associated to a (black 
hole or cosmological) horizon as the (Hawking) temperature associated to the 
radiation emitted from that horizon: a quantum effect as explicitly shown by 
the famous formula Th = hc^ K/2iTkB- As long as we limit to consider station- 
ary Killing horizons all the derivations of this formula are basically equivalent; 
but, going from stationary to dynamical geometries, things make tough mostly 
due to the lack of a time translation symmetry generator. Limiting to spherical 
symmetric geometries, it is possible to show [21 the existence of a vector field 
which, defined through the condition (/C"Gq^)'^ = 0, resembles most of the 
amazing Killing vector properties. We shall term this vector the Kodama vector 
field and refer the interested reader to |10H13j for more accurate discussion. By 
choosing Kodama observers as privileged and using semi-classical techniques, 
we can still establish a correspondence between surface gravity and tempera- 
ture in dynamical, spherical symmetric, black hole or cosmological space-times. 
Semi-classical methods such as the tunneling method or the Hamilton-Jacobi 
method have proven so far to be reliable in all the testable conditions (for a 
recent review on the subject [14] and references therein). However, a compar- 
ison between such methods and standard quantum field theoretic calculations 
in dynamical spaces is still something deserving to be done. With the purpose 
of partially filling this gap, we try to understand how far we can push the cor- 
respondence between surface gravity and temperature in dynamical spherically 
symmetric spaces. We shall do this by searching whether a point-like detector 
will register a quasi-thermal excitation rate of the form 

F - Eexp{-2TrE/KH{t)) 

where Knit) is the dynamical horizon surface gravity to be introduced below, 
and E is the energy of the detector's quantum jump. We shall take this feature 
as a hint that the detector feels the vacuum as a mixed quasi-thermal state 
with an effective temperature parameter Tnit) — K/f(i)/27r. The analysis will 
be done for conformally coupled scalar fields, the best approximation we know 
to massless radiation, but we think the conclusions to be drawn will have general 
validity. Important examples will cover black holes, de Sitter space and generic 
FRW space-times. 

The paper is organized as follows. In Section 2, we briefly review basic 
predictions of the tunneling method. Section 3 contains a discussion on the 
Unruh-DeWitt detector which is introduced in order to confirm through quan- 
tum field theoretic calculations the results of previous tunneling papers. In Sec- 
tion 4 we show how the formalism of the previous section applies to a generic 
static black hole, and finite time effects are investigated. Section 5 is devoted 
to FRW space-times, and to a discussion of the realistic model of a universe 
filled with matter and cosmological constant (fim, S^a), where the main conclu- 
sions are presented as consequence of analytic results supported by numerical 
computations. Concluding remarks follow at the end of the paper. 

We use the metric signature (— , -I-); Greek indices run over to 3 while 
mid-Latin as i,j only over and 1. We use Planck units in which c = h = G = 
kB = 1- 
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2 The Hamilton-Jacobi method 

In 2000 Parikh and Wilczek [15] (see also [TB]) introduced the so-called tun- 
neling approach for investigating Hawking radiation. Here, we shall focus on a 
variant of their method, called Hamilton-Jacobi tunneling method [TT^BD] . This 
method is covariant and can be extended to the dynamical case PDHT^U^ , and 
to the study of decay of massive particles and particle creation by naked sin- 
gularities |22J . In their approach, Parikh and Wilczek made a clever use of the 
Painleve-Gullstrand stationary gauge for four-dimensional Schwarzschild black 
hole 



ds^ = - 1 ]dtl + 2\ drdt^ + dr^ 



r I '' V r ^ 



+ r'^{de'^ +svc?ed(jp-) (2.1) 

which is regular on the trapping horizon ru — 2M 0. This is one of the key 
points since the use of singular gauges, as the Schwarzschild gauge, leads, in 
general, to ambiguities and it is useless in the dynamical case which we are 
interested in. The second merit we can address to Parikh and Wilczek work 
is the treatment of back-reaction on the metric, based on energy conservation. 
In the following, we shall limit to leading term results and neglect the issue 
of back-reaction. However, it may be worth to recall that in the limit where 
the number of emitted quanta is reasonably large back-reaction effects can be 
accounted for by assuming that the mass parameter M is a continuous function 
of time tp. Penrose's diagrams for this more general case have been determined 
too, e.g. in [2511^. 

The Hamilton-Jacobi method is reasonably simple, even though subtleties 
are present (see, for example |1S]). It is based on the computation of the classical 
action / along a trajectory starting slightly behind the trapping horizon but 
ending in the bulk, and the associated WKB approximation (keeping h explicit) 

Amplitude oc e*^ . (2.2) 

For an evaporating black hole such a trajectory would be classically forbidden 
since at the trapping horizon photons are only momentarily at rest, whence 
dr/dtp < inside the horizor^. The related semiclassical emission rate reads 

F cx lAmplitudep cx e^^^ . (2.3) 

with 3 standing for the imaginary part. In the tunneling across the horizon, the 
imaginary part of the classical action / stems from the interpretation of a formal 
horizon divergence, and in evaluating it, one has to make use of Feynman's 
prescription related to a simple pole in integration. This corresponds to the 
choice of suitable boundary conditions in quantum field theory approach. 

It should be mentioned also that in the static/stationary case, there exist 
different interpretations to the Hamilton-Jacobi method (see, for example [2S1), 
none of them, however, can be easily extended to the dynamical case, the one 
we are mainly interested in. 



^ The global event horizon is rg < ru for an evaporating black hole and ro > r[f for an 
accreting black hole, as can be seen from the equation of radial null rays, ro = 1 — ^th/to. 
^ It means the photon must go back in time tp to escape the horizon. 
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We may anticipate that in the WKB approximation of the tunnehng prob- 
abihty, one asymptotically gets a Boltzmann factor, 

Tcxe-^", (2.4) 

where w represents the Kodama energy of the tunneling particle, w = —JC ■ dl. 
It is crucial in our approach that the argument of the exponent be a coordi- 
nate scalar (invariant quantity), since otherwise no physical meaning can be 
addressed to T . In particular, in the Boltzmann factor, /? and the energy oj have 
to be separately coordinate scalars, otherwise again no invariant meaning could 
be given to the quantity P/h. In the static case, we interpret T = h/ P as the 
horizon temperature. The tunneling method for spherically symmetric dynam- 
ical case, has been considered in [H1[ITJ[T^[371[1S] , where the Kodama-Hayward 
invariant formalism has been used. 

In the cosmological case on the other hand, one can still have a trapping 
horizon despite the absence of collapsed matter simply as a result of the expan- 
sion of the universe. As for an excreting black hole, this too is represented by a 
time-like hypersurface. Similarly, an approximate notion of temperature can be 
associated to such horizons based on the existence of a surface gravity and again 
the tunneling method gives a non vanishing amplitude having the Boltzmann 
form. However, we will see that a co-moving monopole detector seems to react 
to the expansion in a different, non "Boltzmannian" way, while reaching ther- 
mal equilibrium (or better, detailed balance conditions) only asymptotically for 
large times. 

As showed in [T2], in a generic FRW space-time, one has 

r - exp(-29/) ~ e"(^'^" , (2.5) 

where, since kh < and ojh > for physical particles, the imaginary part of the 
action is positive definite. This result is invariant since the quantities appearing 
in the imaginary part are manifestly invariant. Furthermore T = — K///27r 
satisfies a First Law. As a consequence, at least in some asymptotic regime and 
for slowly changes in the geometry, one could interpret T — — K///27r as the 
dynamical temperature associated with FRW space-times. In particular, this 
gives naturally a positive temperature for de Sitter space-time, a long debated 
question years ago, usually resolved by changing the sign of the horizon's energy. 
It should be noted that in literature, the dynamical temperature is usually given 
in the form T = H/2-k (exceptions are the papers [211130]) with = A/3 = 
-ffg. Of course this is the expected result for de Sitter space in inflationary 
coordinates, but it ceases to be correct in any other coordinate system. This 
fact seems not so widely known so, for sake of completeness, we shall try to 
show it in detail, de Sitter space in the global patch is described by the metric 

ds^ = -dt^ + a^{t)dnf^^ (2.6) 

with a{t) — cosh(iJoi), '^^^(3) the unit three-sphere and k = Hq is the only 
relevant scale. The Hubble parameter is time dependent, H{t) = Hq tanh(_ffoi), 
and satisfies the identity H{t) = Hq — H^{t) = k/a'^{t). The horizon radius 
turns out to be 
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and Hayward's surface gravity (minus sign due to conventions) is 

^H = - (^H\t) + ]^H{t) + RH{t) = Ho (2.8) 

an invariant quantity, indeed. Hence, we see that the H and k terms have to be 
present in a generic FRW space-time. The important spatially flat case straight- 
forwardly follows, kh = —[H{t) + H{t)/2H{t)]. Note that this is independent on 
position, suggesting that kh really is an intrinsic property of FRW space linked 
to the bulk. The horizon's temperature and the ensuing heating of matter was 
foreseen several years ago in the interesting paper |31j . 

3 The Unruh-DeWitt detector 

We recall that for the decay rate of a massive particle in de Sitter space, the 
exact quantum field theory calculation of Moschella et al. [32] supports the 
WKB semiclassical tunneling result of [2^155] . 

What about the energy scale associated with the horizon tunneling? We 
have shown that the semiclassical WKB method leads to an asymptotic particle 
production rate, involving the Boltzmann factor and related "temperature" 

r~e"^"" ^r = ^ (3.1) 

27r 

where kh is the Hayward invariant surface gravity. For a generic spherically 
symmetric space-time, this result obtained by the Hamilton-Jacobi tunneling 
method seems a very clear prediction, namely an answer to the question: how 
hot is our expanding universe? A possible way to understand this issue us- 
ing quantum field theory in curved space-time is to make use of a "quantum 
thermometer" (basically, a Unruh-DeWitt detector) and evaluate its response 
function, that is, loosely speaking, the number of clicks per unit proper time it 
detects as it is carried around the universe. For a recent review, see [33]. 

In our approach, since we would like to obtain an invariant result, we will 
consider detectors which follow Kodama trajectories in a generic spherically 
symmetric space-time. The problem of back-reaction on these Kodama trajec- 
tories has been investigated in [55] . 

As we will see, the Unruh-DeWitt thermometer gives a clean answer only 
in the stationary case, while for FRW space-time the situation is not so simple, 
since horizon effects are entangled with highly non trivial kinematic effects. For 
general trajectories in flat space-time see the recent paper |36| . An interesting 
analysis has been also put forward by Obadia [ST]- In a recent paper [3^, local 
scaling limit techniques have been used in investigating the Hawking radiation. 

In the following, we review the well know Unruh-DeWitt detector formal- 
ism, adapted to a spherically symmetric conformally flat space-time, namely, 
introducing the conformal time t] by means of dr/ = ^ , the flat FRW space-time 
we are going to deal with reads 

ds'^ = a^{-q){-djf + dy?) , a; = (r;, x) . (3.2) 

For the purpose at hand, it is very convenient to consider a free massless scalar 
field which is conformally coupled to gravity, since, as is well known, the related 
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Wightman function W{x, x') can be computed in an exact way. In fact, one has 

W^(x,x')-^/k(x)/^(x'), (3.3) 

k 

where the mode functions fk{x) satisfy the conformahy invariant equation {TZ 
being the curvature scalar) 

□-?Vk(^)=0. (3.4) 



Making the ansatz 
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Mx) - ^e-^'^- , (3.5) 



one has, for the unknown quantity g{k), 

9"{V) + k'giv) = . (3.6) 
As a consequence, making the choice of the vacuum given by 



one has [5] 



1 1 /47r2 

W{x,x')^ (3.8) 

a{r])a(r]' ) |x — x'p — if] ~ rj' — ie\'^ 

As is usual in distribution theory we shall leave understood the limit as e 0+. 
However, it has been shown by Takagi [35] and Schlicht [301 that this prescription 
is manifestly non-covariant. Since one is dealing with distributions, the limit 
e — > 0+ has to be taken in the weak sense, and it may lead to unphysical results 
with regard to instantaneous proper-time rate in Minkowski space-time. We 
adapt Schlicht 's proposal to our conformally flat case, namely 

1 1 /47r2 

W{x,x')^ , w n r? n . r a, (3-9) 

a[r])a{r]') [[x ~ x' ) — ie[x + x' j]"' 

where an over dot stands for derivative with respect to proper time In the 
flat case, this result has been generalized by Milgrom and Obadia, who made 
use of an analytical proper-time regularization [UlU^]. It should be noted the 
appearance of Minkowski contribution, as a function of the conformal time rj. 

The transition probability per unit proper time of the detector depends on 
the response function per unit proper time which, for radial trajectories, at finite 
time r may be written as |43| 

dF 1 r-^" , e-'^" 
rRe / as 



dr 2Tr^ Jq a{T)a{T — s) 

[x{t) — xir — s) ~ ie{x{T) + xir — s))Y ^ ^ 



where tq is the detector's proper time at which we turn on the detector, and 
E is the energy associated with the excited detector state (we are considering 
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E > 0). Although the covariant ie-prescription is necessary in order to deal with 
the second order pole at s = 0, one may try to avoid the awkward limit e — >■ 0"*" 
by omitting the e-terms but subtracting the leading pole at s = (see [13] for 
details). In fact, the normalization condition 

g^^x'^x'' = [aiT)x{T)f ^-l, (3.11) 

characteristic of time-like four- velocities, has to be imposed. Thus, introducing 
the notation 

(t^{t, s) EE a{T)a{T - s)[x{t) - x{t - s)f , (3.12) 

due to (|3.11|) . for small s, one has 

a^{T, s) = -.s^[l + s^d{T, s)] . (3.13) 

As a consequence, for At ~ t — tq > 0, one can present the detector transition 
probability per unit time in the form 

dF 1 /I 1 \ 

where the "tail" or finite time fluctuating term is given by 

1 cos(i;5) 

The convergence at infinity is assumed, but in all physically interesting cases it 
is ensured. This is not quite the original expression found in |43) but can be 
obtained from it by simple manipulations. This is the main formula which we 
will use in the following. Equation (j3.14l) may be much more convenient to deal 
with than the original expression containing the e-terms, since in the latter the 
limit in distributional sense must also be taken at the end of any computation. 
In the important stationary case in which a{T, s)^ = o'^(s) — cr^(— s), Eq. p.l4p 
simply becomes 

In this case (examples are the static black hole and the FRW dc Sitter space) 
all the finite time dependence is contained in the fiuctuating tail. 



As a result we have the manageable expression (j3.14p . in which the last 
fluctuating tail term incorporates part of the finite-time effects and, as we will 
see, in the case of asymptotically stationary situations controls how fast the 
thermal equilibrium is reached. 



4 Quantum thermometers in static and stationary spaces 

As an application of the formalism previously developed, first we are going 
to revisit the case of the generic static black hole, then we shall consider the 
perturbation of thermal equilibrium due to finite-time effects. 
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4.1 The generic static black hole 

The general metric for a static black hole reads 

ds^ = -V{r)df + + r^dil' , (4.1) 

W{r) 

where, for sake of simplicity, we shall assume W{r) = V{r), with V{r) having 
just simple poles in order to describe what we might call a nice black hole 
(suggested by Hayward). Let th be the (greatest) solution V{r) = 0, the 
general formalism tells us that the horizon is located at r = th] the Kodama 
vector coincides with the usual Killing vector (1,0); and the Hayward surface 
gravity is the Killing surface gravity, namely kh — k = Vj^ /2. We now introduce 
the Kruskal-like gauge associated with this static black hole solution. The first 
step consists in introducing the tortoise coordinate 

^^M^TtI-T- (4-2) 



V{r) 



One has — oo < r* < oo and 



ds^ = V{r*)[-dt^ + {dr*f]+r^{r*)dn%) . (4.3) 



The Kruskal-like coordinates are 



so that 



R = ie"''* cosh(Kt) , T = -e'^''* sinh(Ki) , (4.4) 

K K 



T^ + R^^ ^e""'-* , (4.5) 



and the line element becomes 



ds^ = e-^'^''" V{r*)[-dT^ + dR^]+r^{T,R)dn^ 

= e'^'^'-'\-dT^ +dR^)+r^{T,R)dn^ (4.6) 

where now the coordinates are T and R, r* — r*{T, R), e**^*"*' :— F(r*)e^^'"' . 

The key point to recall here is that in the Kruskal gauge (j4.6p the normal 
metric - the important one for radial trajectories - is conformally related to two 
dimensional Minkoswki space-time. The second observation is that Kodama 
observers are defined by the integral curves associated with the Kodama vector, 
thus the areal radius r(T, R) and r* are constant. As a consequence, the proper 
time along Kodama trajectories reads 

dr^ = Vir*)dt^ = e*('^')(dT2 - dR^) 

= a^{r*){dT^ -dR^) (4.7) 



so that t ~ T/^yV{r*) and 
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The geodesic distance reads 
and one gets, using (j4.8p . 



(r(r) - T(r - s)) V 



+ {R{T)-R{r-s)y 



■ sinh 



(4.9) 



(4.10) 



Since cr^(r, s) = cr^(s) = ct^(— s), we can use equation (j3.16p in the hmit when 
At goes to infinity: 



d7 



1 



1 



sinh X 



(4.11) 



The integral can be evaluated by the theorem of residues and the final result is 

E 



dF_ 
dT 



1 
2^ 



exp 



(^) 



(4.12) 



Since the transition rate exhibits the characteristic Planck distribution, it means 
that the Unruh-DeWitt thermometer in the generic spherically symmetric black 
hole space-time detects a quantum system in thermal equilibrium at the local 
temperature 



T 



(4.13) 



With regard to the factor y/V* — ^/~goo, recall Tolman's theorem which states 
that, for a gravitational system at thermal equilibrium, Ty/—goQ = constant. 
For asymptotically flat black hole space-times, one obtains the "intrinsic" con- 
stant temperature of the Hawking effect, i.e. 



T„ = — 



2n An 



(4.14) 



We would like to point out that this is a quite general result, valid for a large 
class of nice black holes, as for example Reissner-Nordstrom and Schwarzschild- 
AdS black holes. On the other hand, the Schwarzschild-dS black hole cannot 
be included, due to the presence of two horizons. However, as an important 
particular case, we may consider the static patch of de Sitter space, with a 
metric defined by 



v{r) = 1 - ny , 



A 



(4.15) 



The unique horizon is located at rn — Hq^ and the Gibbons-Hawking temper- 
ature is [35 Th = Hq/2'k. In the next Section, we will present a derivation of 
this well known result in another gauge. 
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We conclude this subsection making some remarks on de Sitter and anti-de 
Sitter black holes. First, we observe that in a static space-time, namely the one 
corresponding to a nice black hole, the Killing-Kodama observers with r — K 
constant have an invariant acceleration 

A- ^ g,.A^A^ = ^ , (4.16) 

where — u^Vi/u'*, being the observer's four-velocity, that is the (normal- 
ized) tangent vector to the integral curves of the Kodama vector field. In the 
case of de Sitter black hole, one has 



^0 

As a result. 



and the de Sitter local temperature felt by the Unruh detector 

Hp 1 
2^ VI - H^K^ 



T<is = ^^=^— (4.19) 



can be re-written as 



Tds = ^ V^' + ^0 = + Thu ■ (4-20) 



A similar result was also obtained for AdS in 46 , and it reads 

1 



TAdS = ^^/A^~Hi. (4.21) 

We would like to show that it is a particular case of our general formula (|4.13l) . 
In fact, it is suHicient to apply it to the four-dimensional topological black hole 
with hyperbolic horizon manifold found in [I5H5T| . which is a nice black hole 
with 

Vir) = -l-- + Hy , (4.22) 
r 

where C is a constant of integration related to the black hole mass. The space- 
time is a solution of Einstein equation with negative cosmological constant A = 
—3Hq, which is asymptotically Anti-de Sitter. When the constant of integration 
goes to zero, one has still a black hole solution, and calculation similar to the 
one valid for de Sitter space-time gives 



which is Deser et al. result [IB]. Thus, for spherically symmetric space-times 
with constant curvature one has that the local temperature felt by the Kodama- 
Unruh-DeWitt detector can be written as 



T^Jn+^nH, (4.24) 
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where Tjj is the Unruh temperature associated with the acceleration of the 
Kodama observer, Tqh is the Gibbons-Hawking temperature and a = 1 for 
the de Sitter space-time, a ~ Q for Minkowski space-time (this is the original 
Unruh effect) and a = — 1 for the "massless" AdS topological black hole. This 
formula may help to understand better the relation between the Unruh-like 
effects and the genuine presence of a thermal bath and shows that, in general, the 
Kodama-Unruh detector gives an intricate relation between Killing-Hayward 
temperature and other invariant temperatures such as the Unruh's one. Note 
that T in Eq. (j4.13p is greater than Tu = A{r)/2TT for r > rn, where A is the 
local acceleration of an observer following a Killing trajectory in the black hole 
space-time, a fact that has been interpreted as a violation of the equivalence 
principle [17] . We prefer to interpret this effect as due to the additional presence 
of the Hawking radiation over the vacuum thermal Unruh's noise. 



4.2 Finite-time effects in stationary space-times 

We now present a brief discussion of finite-time effects which will be relevant 
to the following discussion on asymptotic behaviour: how is the thermal distri- 
bution of the response function reached in the limit of very large times? In the 
case of non inertial particle detector in Minkowski space-time, see |52] , and for 
de Sitter FRW space see [55] . 

To answer this, we consider the finite time contribution due to the fluctuating 
tail (the Jr term in Eq. p.lGI) ') for de Sitter or black hole cases compared to the 
thermal value given by the time-independent part. A direct calculation of the 
tail (I3.15|) using Eq. (|4.10|) for black holes (in particular Eq. (|4.15p for dS) and 
the fact that 

oo 

csch2(a;) =4^ne-2"" (4.25) 



gives 



n=l 



• cos{Es) 
B7r2 sinh^ {^) 

2^ n2-H£;2/47r2T^ ^ 

X rMk„cos(i;Ar) - sin(£;AT) ) (4.26) 



where ni is the rescaled surface gravity and Th = ki/2t: = k/2t:\'V the local 
Hawking temperature. We recall that k = for de Sitter space and k — 
V'{rH)/2 for the black hole: these quantities in fact determine the characteristic 
time-scales the thermalization time has to be compared to. 

We consider as before the peculiar Kodama observer for which V — 1, so 
that Th — k/2tt. As a general feature, the fluctuating tail term drops out 
exponentially for large At, that is for long proper time intervals in which the 
detector stays on. In order to analyze the approaching to an equilibrium con- 
dition of the response function, we consider the ratio between the finite-time 
expression - the sum of F and the tail Jr - and F alone, with the agreement 
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that equilibrium is attained whenever 

Re, - - 0(1) . 

r 

Looking at (I4.26P one easily sees that the equilibrium value, Req = 1, is reached 
sooner ii E/ k ^ 1. To be more precise, irrespective of the absolute value of n, 
a detector that is switched on for a time much shorter than the characteristic 
time-scale At <C k^^, detects a thermal bath only with particles whose energies 
are E <C k; on the other hand, a thermal equilibrium for particles with energies 
E ^ K is registered only if the detector lifetime is At 3> which is the age 
of the universe. The Hubble scale corresponds to an extremely small energy 
scale of order lO^^^Gev, therefore E ^ k is the physical region. 

It easy to see that if the factor V < 1, the thermalization time decreases for 
every energy scale. 



5 The FRW and asymptotically de Sitter space-times 

To apply the Unruh-DeWitt detector formalism to cosmology we consider a 
generic FRW spatially fiat space-time. This case has been investigated also 
in [37] (see also [M]). Recall that here the areal radius is i? = ra{t). Thus, for 
the Kodama observer, one has 

with constant K. For a radial trajectory, the proper time in FRW is 

dT^ = a{r]){d7]^ - dr^) . (5.2) 

As a function of the proper time, the conformal time along a Kodama trajectory 
is 

V{r) ^ - [ dr 



= - [ dr , (5.3) 

H{t) being the Hubble parameter as a function of proper time. In general, we 
may use Eq. p.l4p in which, for radial Kodama observer, one has 



x{t) = (77(r),r(T),0,0) 

1 K \ 
; dr,— -,0,0 . (5.4) 

As a warm up, we first revisit the well known example of FRW space is the 
stationary flat de Sitter expanding (contracting) space-time, which in the FRW 
context is defined by considering a{t) — e^"*. Thus, 

ds^ ^ -dt^ + e^"°\dr^ +r^dn^) . (5.5) 



5 The FRW and asymptotically de Sitter space-times 



13 



Here H{t) = Hq is constant as well as ^ = Vb = 1 — HqK^. For Kodama 
observers 



air) 



and 



r]{T) = e , rir) — K e ^ 

Ho 



thus, the geodesic distance is 



Hq 



Hps 
2V% 



(5.6) 



(5.7) 



(5.8) 



This result is formally equal to the one obtained for the generic static black hole 
(|4.10p . Since again (T^(r, s) = cr^(s) = (t^(— s), we may use (I3.16|) and obtain, 
for £' > and in the infinite time limit 







dF _ i/( 



dx ( 



Again, we arrive at 



dFds ^ 
dr 27r 



E 



1 



1 



sinh X 



exp 



V Ha 



1 



(5.9) 



(5.10) 



which shows again that the Unruh-DeWitt thermometer in the FRW de Sit- 
ter space detects a quantum system in thermal equilibrium at a temperature 
T = Hp/2-Ky/VQ. Here, the Tolman factor takes the form a Lorentz 7-factor, 
which represents the Unruh acceleration part. In fact, we recall that the four- 
acceleration of a Kodama observer in a FRW space-time has the expression 



A' 



H{t) + {1- H^{t)K^)H^t) 



As a result, for dS space in a time dependent spatially flat patch we have 



(5.11) 



A^ 



1 - K^Hl ' 



(5.12) 



showing that 



Ho 



Hi + 



in agreement with the dS static calculation. When X = 0, that is when the 
detector is co-moving, one has Vq = 1 and the classical Gibbons-Hawking result 
Tds — H0/2-K is recovered. 



Let us come to consider the more realistic scenario of a truly dynamical 
space-time of cosmological interest. From previous considerations, our basic 
formulas for the transition rate of the detector p.l4p are manageable - in the 
sense that we can extract quantitative information - only in the few highly 
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symmetrical circumstances mentioned in Section 4. As it will be clear at the 
end of this section, any departure from those models is responsible for significant 
difficulties. For instance, let us take on the case of homogeneous, spatially flat, 
universe dominated by cold matter and cosmological constant. The scale factor 
is (e.g. [55]) 



i(<) ^aosinh^/^ Qv^^^o 



t] (5.13) 



where oq = {flm/^Ay^^ and flm + f^A = 1; -f^o — \/87rpcr/3 and VIa rep- 
resents the relative density of matter (if A = m) or cosmological constant (if 
A — A). Setting h = \/fl\HQ for simplicity, its current value is of order 
/i w 2 X lO^^^sec""'^. Upon integration the conformal time becomes 




where 2Fi{a,b,c; z) is a hypergeometric function and the constant has been 
opportunely chosen so that at the Big Bang ri{t = 0) = 0. The detector's proper 
time is related to the cosmic time through a manageable expression only if we 
limit ourselves to consider co-moving detectors: T(i) - tq = / dt^/l - H'^(t)K'^ 
so that for = 0, AT(t) = i. At being the proper time interval during which 
the detector is turned on. Unlike the stationary cases analyzed previously, this 
model presents a Big Bang singularity at the origin of the time coordinate, so 
that the detector must be switched on at some tq > 0. In particular, the Big 
Bang prevents taking the limit as tq — > —00. By the same reason, the scale 
factor (|5.13p is defined only for positive values of the argument: a new feature 
with respect to what we have seen in the previous Sections. As a consequence, 
a{t — s) is defined as in (j5.13p only for i — s > and trivially continued outside 
the interval in order to make well defined the transition rate p.l4[) . 

Leaving the technical details to the Appendix A, we obtain the following 
response function 

Ft — FdS + JdS,T + 



27r' 

n=l fe=l 



^hkAr 



hk cos(£'At) + E smiEAr) ) - hk 

(5.15) 



h^k^ + E^ 

in which F^s is the De Sitter r independent contribution given by Eq. (j5.10l) but 
with effective Hubble constant h = \/TIaHo and Jds,T is the related tail given 
by Eq. (j4.26p . The numerical coefficients g{n,k) can in principle be computed 
but enter in a tail which decays exponentially fast in the switching time and 
which also contain oscillating terms. 

We may take the limit At — > 00 and observe that every t dependent term 
of this expression goes to zero. This is the main result of our paper. To 
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summarize, we may say that the detector cUcks close to a de Sitter response 
and reaches thermahzation (possibly, through decaying oscillations) as Ar is 
sufficiently large. In fact, as far as the regime hAr ^ 1 is concerned, de Sitter 
space-time is recovered. We may think of this as describing a de Sitter thermal 
noise continuously perturbed by the expansion (or contraction) of the universe. 
In particular, insofar as we can speak of temperature, in this large-time regime 
the detector registers the de Sitter temperature h/2n, equal to the large-time 
limit of the horizon temperature parameter given by the surface gravity, which 
in the present case has the exact but slow long-time evolution 

Th = — [coth(3/it/2) - 3/4sech(3/it/2)csch(3/it/2)] 
27r 

It is worth noting that, while in the stationary phase hAr ^ 1 the limiting 
result is consistent with the limiting value of the surface gravity, in the non- 
stationary regime it seems less trivial to compare the results of the two methods, 
because it has not been possible to extract a temperature parameter from the 
transition rate of the detector, but asymptotically. 

6 Conclusions 

In this paper, with the aim to better understand the temperature-versus-surface 
gravity paradigm, the asymptotic results obtained by semiclassical method in 
previous papers have been tested with more reliable quantum field theory tech- 
niques as the Unruh-DeWitt detector analysis. For black holes and pure de 
Sitter space the two analysis are mutually consistent and even predict the de- 
pendence of the temperature on position or acceleration. Moreover, the analysis 
of the oscillating tail has been extended to stationary black holes. 

For cosmology and away from de Sitter space the thermal interpretation, 
strictly speaking, is lost but the detector still gets excited by the expansion of 
the universe. By accepting the surface gravity versus temperature paradigm we 
would expect a quasi-thermal excitation rate of the form 

F ^ Eex^pi-E/Tnit)) 

Tnit) being given by our last expression above. That is, although in a generic 
FRW space-time the thermal interpretation breaks down in most of the cases 
because of the time-dependence of the background, still this time dependence of 
the transition rate could be expected to mainly reside in an effective temperature 
parameter. But using a comoving detector this is not what we have found. For 
instance, in the Einstein-de Sitter regime there seems to be excitations of non- 
thermal type and we showed that the scale factor of the flat ACDM-cosmology 
has no other temperatures in action than the de Sitter one. 

It remains to see whether there is any non trivial quasi-thermal effect on 
accelerated, or more general, Kodama trajectories. In the affirmative case, 
that would mean that the horizon surface gravity and temperature should be 
associated more likely to vacuum correlations than to particle creation and 
forces, in our view, a different interpretation of the tunneling picture. In this 
respect, the classical Parker's papers on particle creation [SHIIST] are certainly 
relevant. One possibility is that the horizon surface gravity could represent 
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an intrinsic property of the horizon itself, leading to some kind of holographic 
description, while the detector in the bulk simply clicks because it is embedded 
in a changing geometry. In fact, we would expect the clicks in almost any 
changing geometry, even for those lacking a trapping horizon. 
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APPENDIX A: Response function for ACDM model 

In order to obtain (|5.15|) . we define the variable x = cxp (— |/iAr) and expand 
the inverse cr'^{x, s), given by the scale factor (j5.13p . around x = {i.e. At — > 
c»). We obtain a reasonably simple expansion in even powers of x given by 

oo / 3n— 1 \ 

On the right hand side, the first term is the constant term of the expansion and 
happens to be the pure de Sitter contribution, i.e. 

als{s) = -Ah-'sin\.'{^^ (A-2) 

with the effective Hubble constant h — ^/^IaHq. Numerical hints given by 
the coefficients of the expansion up to the 10*'' order in x, allow us to make a 
conjecture that the g{n, fc)'s in the second term have a mean decreasing behavior 
and are bounded in the interval (0,1), but the main point is that the series 
in ea. (jA-l[) is absolutely convergent with a finite radius of convergence which 
includes any i > 0, namely the entire range of integration. 

Hence, integrating term by term the expression (jA-ip . and making use of 
(|3TT4)) . for finite Ar one has ea. ((5?T5t . 

In the At — cxi limit, we can focus on the leading exponentials contained 
in the last term of this expression: these are the k = (3n — 1) terms, which are 
all dominated by a common factor exp (—hkAr). All the other terms are even 
more damped, so the convergence to zero is evident. 
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